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Abstract. We investigate the appearance of it lapses in the transmission phase 9 
of a two-level quantum dot with Coulomb interaction U. Using the numerical and 
functional renormalization group methods we study the entire parameter space for 
spin-polarized as well as spin-degenerate dots, modeled by spinless or spinful electrons, 
respectively. We investigate the effect of finite temperatures T. For small T and 
sufficiently small single-particle spacings 6 of the dot levels we find it phase lapses 
between two transmission peaks in an overwhelming part of the parameter space of 
the level-lead couplings. For large S the appearance or not of a phase lapse between 
resonances depends on the relative sign of the level-lead couplings in analogy to the 
U = case. We show that this generic scenario is the same for spin-polarized and 
spin-degenerate dots. We emphasize that in contrast to dots with more levels, for 
a two-level dot with small S and generic dot-lead couplings (that is up to cases with 
special symmetry) the "universal" phase lapse behavior is already established at U = 0. 
The most important effect of the Coulomb interaction is to increase the separation of 
the transmission resonances. The relation of the appearance of phase lapses to the 
inversion of the population of the dot levels is discussed. For the spin-polarized case 
and low temperatures we compare our results to recent mean-field studies. For small 
6 correlations are found to strongly alter the mean-field picture. 
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1. Introduction 

The local Coulomb interaction U > of electrons occupying quantum dots leads to a 
variety of effects. Many of them can conveniently be studied in transport through the 
dot within the linear regime. Theoretically as well as experimentally well-investigated 
examples are the Coulomb blockade (CB) peaks of the transmission (conductance) [1] 
as well as the plateaus of width U of the transmission (conductance) induced by the 
Kondo effect [2j. Additional features of interacting multi- level dots that have recently 
attracted considerable theoretical attention are the population inversions of the dot 
levels [31 HI El E] , the phase lapses of the transmission phase 9 or, equivalently, the zeros 
of the transmission amplitude t (transmission zeros) [3J, [Tl [8], [9j [101 [11] and additional 
correlation- induced resonances of \t\ [12]. They appear in certain parts of the parameter 
space when the level occupancies and the transmission amplitude are investigated as 
functions of the level positions, which can be tuned via a nearby plunger gate voltage. 
Such effects were hitherto mostly studied in a minimal model involving only two levels. 
A very important step towards a unified understanding of population inversions, phase 
lapses and correlation-induced resonancess in spin-polarized two-level dots was recently 
taken by a multi-stage mapping of the problem on a generalized Kondo model and a 
subsequent renormalization group and Bethe ansatz analysis of the effective Hamiltonian 

Theoretical studies of phase lapses (transmission zeros) are of primary interest in 
connection with a series of linear response transmission measurements by the Weizman 
group [HI [15], [16] on Aharonov-Bohm rings containing a quantum dot in one arm. Under 
suitable conditions both the phase 9 and magnitude \t\ of the transmission amplitude 
t = \t\e l9 of the dot can be extracted from the Aharonov-Bohm oscillations of the current 
through the ring [T7]. When this is done as function of a plunger gate voltage V g that 
linearly shifts the dot's single-particle energy levels downward, Ej — — V g (j — 1, 2, . . . 
is a level index), a series of well-separated CB peaks of rather similar width and height 
was observed in |t(V^)|, across which 9(V g ) continuously increased by ir, as expected 
for Breit-Wigner-like resonances. In each CB valley between any two successive peaks, 
9 always jumped sharply downward by 7r. This phase lapse behavior was found to be 
"universal" , occurring in a large succession of valleys for every man?/-electron dot studied 
in Refs. [HI [151 [IS]- This universality is puzzling, since naively the behavior of 9(V g ) is 
expected to be "mesoscopic" , i.e. to show a phase lapse in some CB valleys and none 
in others, depending on the dot's shape, the parity of its orbital wavef unctions, etc. 
Only recently [16], also the /eiy-electron regime was probed experimentally: as V g was 
increased to successively fill up the dot with electrons, starting from electron number 
N e = 0, 9{V g ) was observed to behave mesoscopically in the Jew-electron regime, whereas 
the above-mentioned universal phase lapse behavior emerged only in the many-electron 
regime (N c > 15). 

It was suggested in Ref. [16] that a generic difference between the few- and many- 
electron dots may be that for the latter, transport might simultaneously occur through 
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several partially filled single-particle levels in parallel. A possible reason could be that 
the mean (non-interacting) level spacing S of the topmost filled levels decreases as the 
number of electrons increases, while the charging energy U still implies well separated 
transmission resonances [IB] . This scenario forms the basis of a recent systematic 
study by us of the interplay of level spacing, level width and charging energy on the 
phase lapses for up to four interacting levels and spin-polarized electrons [19]. We 
showed that the universal phase lapse and transmission zero behavior appearing at 
small 5 can be understood as resulting from a Fano-type interference effect [20] involving 
transport through two or more effective dot levels, whose position and width have been 
renormalized by the Coulomb interaction and coupling to the leads. The importance of 
several overlapping levels for phase lapses had earlier been pointed out by Silvestrov and 
Imry in a rather specific model of a single wide and several narrow levels with strong 
interaction [3]. 

Here we supplement our earlier study [19] by discussing the relation between phase 
lapses and population inversions and by investigating the role of finite temperatures 
T > as well as spin, focusing on N = 2 levels. When spin is included, the Kondo 
effect plays a role for an odd average occupation of the dot, but the phase lapse scenario 
is unaffected by this. Experimentally the behavior of the phase in the presence of the 
Kondo effect was investigated in Refs. [21] and [22]. As in Ref. [19] we are concerned 
with the generic behavior and thus investigate the entire parameter space, going beyond 
subspaces of higher symmetry [such as left-right (1-r) or 1-2 symmetry of the coupling 
between the left and right leads and the two levels] . For low temperatures and sufficiently 
small single-particle spacings S of the dot levels, we find 7r phase lapses between two 
transmission peaks in an overwhelmingly large part of the parameter space of the level- 
lead couplings. We point out that the two level case is special compared to models with 
N > 2, as for generic level-lead couplings a transmission zero and phase lapse occurs 
between the two transmission peaks even at U=0. The effect of the interaction is merely 
to increase the separation of the transmission peaks. For large S the appearance or not 
of phase lapses between transmission peaks depends on the relative sign of the level-lead 
couplings in analogy to the noninteracting case [9]. 

For spin-polarized dots we in addition compare our T = results with the ones 
of recent mean- field studies pUi EH]. In these works level- lead couplings beyond the 
subspaces with increased symmetry were studied, and a remarkably more complex 
behavior was found once the symmetries were broken. The importance of considering 
such generic parameter sets was independently pointed out in Ref. [12J. We here 
elucidate how the phase lapse behavior is affected by quantum fluctuations, which are 
expected to be strong in low- dimensional systems. We find that upon including quantum 
fluctuations, the part of parameter space exhibiting universal 7r phase lapses between 
well separated CB peaks becomes larger than suggested by the mean-field study. In 
particular, we do not recover certain peculiar features of the mean-field results of Refs. 
[TO] and [IT] , namely the occurrence, in certain regimes of parameter space, of a phase 
lapse of less than 7r (instead of precisely 7r), accompanied by the disappearance of the 
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corresponding transmission zero PHI E]- These features thus turn out to be artifacts 
of the mean-field approximation, which misses the rather simple scenario for the phase 
lapse behavior of a two-level dot at small 5: for generic level-lead couplings a phase 
lapse and transmission zero between two transmission peaks is already present at U = 0; 
increasing the Coulomb interaction the peaks become well separated while the phase 
lapse and transmission zero remain in the valley between them. 

In the model of a single wide and several narrow levels [3] a relation between phase 
lapses and population inversions was discussed. Therefore, in phase lapse studies quite 
often also the level occupancies rij, j = 1,2, are investigated. We emphasize that the 
generic appearance of a phase lapse and transmission zero even at U — renders the 
two-level model unsuitable for establishing a general relation between phase lapses and 
population inversions, as the latter only appear at sufficiently large U. Furthermore, 
we show that discontinuities of the function of V g are an artifact of the mean- 

field solution (see Refs. [10] and [6]). Within our approaches discontinuities are only 
found for 1-r symmetric level-lead couplings with a relative plus sign of the underlying 
hopping matrix elements and degenerate levels, a case which was earlier identified as 
being nongeneric [121 HB], because the transmission shows only a single peak. 

This paper is organized as follows. In Sec. [2] we introduce our model for the 
spin-polarized and spin-degenerate two- level dot. We discuss the relation between the 
measured magnetic flux dependence of the interferometer's linear conductance and the 
magnitude and phase of the dot's transmission amplitude. The latter can be computed 
from the one-particle Green function of the dot. We present a brief account of our 
techniques to obtain the latter, the numerical (NRG) [23J and functional renormalization 
group (fRG) methods. For an introduction to the use of the fRG to quantum dots see 
Refs. [23] and [25J. In contrast to the implementation of the NRG used in Ref. [19J, 
where we were restricted to 1-r symmetric dots at T = 0, we have now adopted the full 
density matrix (FDM) NRG method of [26], which enables us to investigate dots with 
arbitrary level-lead state overlap matrix elements t l - (with I = L,R) as well as to study 
finite temperatures. In Sees. [3}l5]we present our results of the V g dependence of \t\ and 9. 
First we briefly discuss the noninter acting two-level dot with generic level-lead couplings 
and point out that the phase lapse scenario differs from the one for more than two 
levels. We then investigate interacting, spin-polarized dots, study the relation between 
phase lapses and population inversions and compare to the mean-field results for the 
phase lapses. The issue of continuous versus discontinuous V g dependence of the level 
occupancies ni and n 2 is commented on. Next we study the role of finite temperatures. 
Finally, we consider spin-degenerate levels at small T which implies the appearance of 
Kondo physics at odd average dot filling. Using NRG and fRG we show that the spin 
does not alter the universal phase lapse scenario. Our findings are summarized in Sec. 

El 
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2. The model and methods 

In this section we introduce our model for the two- level dot. We argue that it is the 
energy dependent (effective) transmission amplitude t(u) which one has to compute if 
one is interested in comparing to the measurements of Refs. [141 [T5l [T6] of the magnitude 
of the transmission amplitude and its phase. The amplitude t(u) can be determined 
from the matrix elements of the dot's interacting one-particle Green function. We 
furthermore discuss aspects of the NRG and the fRG specific to our problem. 

2.1. Two-level setup and transmission amplitude 
Our Hamiltonian consists of three parts 

H = -fflead + -Hdot + -^lead-dot • (1) 

The two semi-infinite leads are modeled as noninteracting one- dimensional tight-binding 
chains and for simplicity are assumed to be equal 

oo 

#lead = -T Yl J2^m,^,l Cm + 1 ^ 1 + H ' C -) ■ ( 2 ) 

l=L,R a m=0 

The hopping strength in the leads is r. We use standard second quantized notation 
with I = L,R indicating the left and right leads, where the quantum numbers m and a 
label Wannier states and spin, respectively. The dot is described by 

= J2J2 £ j4a> +\ U J2Y1 (4a> - 1) (4,<Av - '( 3 ) 

a 5=1,2 o,a' 

where the term with j = j' and a = a' is excluded from the sum in the interacting 
part. We define £1/2 = =F^/2 — V g . In experimental systems the inter- and intra-level 
Coulomb repulsion can be expected to be comparable in size and to avoid a proliferation 
of parameters we assumed them to be equal. This assumption is not essential; by 
relaxing it we have checked that our results are robust against inter-level variations of 
the interaction strengths. Finally, the coupling between dot and lead states is given by 

^lead-dot = - E 2 (4 C !>,aA> + R ^ ^ 

l=L,R a 3=1,2 

with real overlap matrix elements tj. For later purposes we define s = sign (if tft% tf) 

and 7 = {rf,rf,rf,rf}/r. 

For simplicity, part of our studies will be performed on a model of spinless electrons, 
for which the spin index will be dropped. The resulting model may be regarded as a 
spin-polarized version of the spinful model obtained if the latter is put in a very large 
magnetic field. 

The experimental two-path interferometer is characterized by wide base regions 
and very narrow point contacts towards the emitter and collector. Electrons can thus 
only pass once through either the upper or lower arm (one of them containing the 
dot) before reaching the collector. It is furthermore reasonable to assume that the 
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voltage drops at the point contacts and the dot remains in equilibrium [TTJ . Under these 
conditions the magnetic flux <p dependent part of the conductance of the Aharonov-Bohm 
interferometer Gab is given by [I7J [27J fT6] 

/oo 
dujf'{uS) \t dot (uj)\ |t rcf I COS [27T0/0O + #ref + #dotM] ( 5 ) 
-oo 

with the flux quantum <fio, the derivative /' of the Fermi function as well as the 
transmission amplitudes tdot = |£dot|e l£,dot of the arm containing the dot and t re f = 
\trei \e ldvci of the reference arm. The transmission tdot is the product of the transmission 
t through the dot and the transmission t Test through the rest of the interferometer arm 
containing the dot. It is reasonable to assume that t res t as well as t re f are only weakly 
energy and gate voltage dependent and thus the V^-dependence of the Aharonov-Bohm 
oscillations of the measured linear conductance of the interferometer is dominated by the 
V g dependence of the magnitude and phase of the transmission amplitude through the 
dot. As usual [9], we compute the energy-averaged transmission phase 9 and magnitude 
\t\ of the dot for a fixed spin direction as the phase and absolute value of 

/oo 
dLjf(u)t(u) , (6) 
-oo 

where the frequency integral represents an energy average, weighted by the derivative 
of the Fermi function, as appropriate for a finite electron temperature in the leads. In 
the limit T — > 0, — /' reduces to a 5-function and t(V g ) is equal to £(//). We here take 
the chemical potential \i = 0. 

Using scattering theory t(u) (for fixed spin direction) can be related to the spin- 
independent matrix elements (in the j = 1,2 indices of the Wannier states) of the dot's 
one-particle retarded Green function Q 

i(u) = 2 (Jrtr?g ltl (u + io) + y^rf q x ^ + «o) 

+ S yT^£ 2il (a; + iO) + S yTtrf &, 2 (a; + iO)J , (7) 

with (after taking the wide band limit; see below) 

r; = vr|t;.| 2 A cad(0)>0, (8) 

where pi ea d(^) denotes the local density of states at the end of each semi- infinite lead. 
Without loss of generality we have assumed that t\ > 0, t\ > and t% = s|tf I- The spin- 
independent dot occupancies rij (per spin direction), that we will also investigate, follow 
from the Green function Qjj by integrating over frequency (or can be computed directly 
when using NRG). For 1-r symmetry of the level- lead couplings and at temperature 
T = 0, t(0) can also be expressed in terms of the spin independent occupancies [E] 

t(0) = sin {[n e - n ]n)e i{n * +n ° )/7T , (9) 

where n e = n\ + n 2 , n Q = for s = + and n e = ni, n = n 2 for s = — , respectively. 
Here we will compute Q in two ways, using both a truncated, that is approximate, fRG 
scheme, and a numerically exact method, the NRG. 
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2.2. The fRG approach 

The truncated fRG is an approximation scheme to obtain the self-energy E (and thus the 
one-particle Green function) and higher order vertex functions for many-body problems 
[28j [291 SO] • As a first step in the application of this approach to quantum dots one 
integrates out the noninteracting leads within the functional integral representation of 
our many-body problem [31]. The leads provide a frequency dependent one-particle 
potential on the dot levels. On the imaginary frequency axis it is given by 

v£$, ta ,(iu) = ]Tt;4g lcad M 5 a>a , , (io) 

I 

where gi ea d(^0 denotes the spin-independent Green function of the isolated semi- infinite 
leads taken at the last lattice site 

As we are not interested in band effects we take the wide band limit. The potential 
then reduces to 

i 

After this step, instead of dealing with an infinite system we only have to consider the 
dot of two interacting levels. 

In the computation of the interacting one-particle Green function projected onto 
the dot system the sum of the dot Hamiltonian with U = and V]*^-, a ,{iuS) can be 
interpreted as a frequency dependent "single-particle Hamiltonian" and in the following 
will be denoted by ho(iu). For the spin-polarized case it is a 2 x 2 matrix in the quantum 
number j = 1,2. Including spin, because of the additional quantum number a =f, |, 
ho(iu) is a 4 x 4 matrix which is block-diagonal in a (spin conservation). As we are here 
not interested in the role of a magnetic field lifting the spin-degeneracy of each level 
the a =| and a —[ blocks are equivalent. The resolvent Go{z) — [z — ho(z)] 1 obtained 
from ho(z) is equivalent to the noninteracting propagator of our two-level many-body 
problem projected on the dot levels. In the generating functional of the one-particle 
irreducible vertex functions we replace Qo(iu) by 

g£(iu) = e(M - A)£oM = e(M - a) [iu - froMr 1 (is) 

with A being an infrared cutoff running from oo down to 0. Taking the derivative 
with respect to A one can derive an exact, infinite hierarchy of coupled differential 
equations for vertex functions, such as the self-energy and the one-particle irreducible 
two-particle interaction. In particular, the flow of the self-energy S A (one-particle 
vertex) is determined by S A and the two-particle vertex 7 A , while the flow of 7 A 
is determined by S A , 7 A , and the flowing three-particle vertex. The latter could be 
computed from a flow equation involving the four-particle vertex, and so on. At the 
end of the fRG flow S A=0 is the self-energy S of the original, cutoff- free problem we are 
interested in [281 129] from which the Green function Q can be computed using the Dyson 
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equation. A detailed derivation of the fRG flow equations for a general quantum many- 
body problem that only requires a basic knowledge of the functional integral approach 
to many-particle physics [31] and the application of the method for a simple toy problem 
is presented in Ref. [30]. For an overview of the application to quantum dots see Refs. 
[21] and [25]. 

We here truncate the infinite hierarchy of flow equations by only keeping the 
self-energy and the frequency-independent part of the two-particle vertex. Higher 
order terms can be neglected if the bare two-particle interaction is not too large. 
By comparison to NRG data this approximation scheme was earlier shown to provide 
excellent results for a variety of dot systems [121 EH 03] • For further comparison see Fig. [7] 
below. The present scheme leads to a frequency-independent self-energy (see below). As 
finite frequency effects (inelastic processes) become important at temperatures T > 0, 
but these are not accurately treated by the level of approximation used here, in the 
present paper we shall show fRG results only for T = 0. The truncation leads to the 
coupled differential flow equations 

3 = 4EE eluJ0+ <vM ii<,vAi ( 14 ) 



and 



OA 2tt 

uj=±A LI 



-^lk',l';k,l = ^ Y Y Y {^t,m'( iuJ )Sn,n'(- iuJ hk',V;m,n% 



A 

m'n':k,l 



=±A m,m' n,n' 

+ < 5m,m'(^ U ) Qn,n'i^ u ) [~~ lk',n';k,l lm',l';n,l + lv,n';k,m 7m',fe';n,l] } 

where k, I, etc. are multi-indices representing the quantum numbers j, a and 

Q\iu)= [£ -1 (W) -£ A ] -1 . (16) 

In the model with spin-degenerate levels each index k, I etc. can take four different values 
j = 1, 2 and a =f, J, which gives 16 equations for E A and 256 for the two-particle vertex. 
For a spin-polarized two-level dot the multi-indices take two values and one obtains 
4 equations for S A and 16 for the two-particle vertex. The number of independent 
equations can be significantly reduced (see below) taking into account the antisymmetry 
of the two-particle vertex and the spin symmetry (for spin-degenerate levels) both being 
preserved by Eqs. (JHJ) and (TI5]) . The initial conditions at A = A — » oo are given by 
Sj", = while 7 1 /° 2 /.i 2 1S given by the bare antisymmetrized two-body interaction. 
In the spin-polarized case the only nonzero components of the two-particle vertex at 
A = A — > 00 are 

7i Ao 2; i, 2 = 7 2 Ao i ;2 ,i = U and 7 i A ° 2;2 ,i = 7 2 A °i ; i, 2 = -U . (17) 
In the model including spin the initial conditions take the form 



7lf,lJ.;lf,lJ. _ U ) 7lf,2T;lT,2T _ ' 7lt,2J.;lt,2| 
72T,2J.;2T,2J. U ' 'U,2|;U,2| U > '1|, 2 T;1|,2T 



All other components which do not arise out of these by permutations (7 A 2-i' 2' = 7 A °2 
and 7^2-1' 2' = — 7^2-2' 1') are zero - The self-energy matrix and thus the one-particle 



,2';1,2 
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Green function is completely independent of the spin direction and in the following we 
suppress the spin indices. 

As already mentioned the present approximation leads to a frequency-independent 
self-energy. This allows for a simple single-particle interpretation of its matrix elements. 
The sum of the S^- and the bare level position correspond to the flowing effective 
level positions, e A = Sj + S^-, while t A = — E^ 2 = — is a hopping between the 
levels 1 and 2 generated in the fRG flow. The fRG formalism then reduces to a set of 
coupled differential flow equations for e A , t A and a few (one in the spin-polarized case 
and seven for spin-degenerate levels) independent components of the two-particle vertex. 
These flow equations can easily be integrated numerically using standard routines. It 
is important to note that although we start out with intra- and inter-level Coulomb 
interactions of equal strengths they generically become different during the fRG flow 
(because of the different 1^). Furthermore, additional interaction terms which are 
initially zero will be generated in the flow. The set of equations significantly simplifies 
if the flow of the vertex is neglected while the results remain qualitatively the same. 
Within this additional approximation and for a spin-polarized dot the flow equations 
for e A and t A are explicitly given in Ref. [12]. In certain limiting cases it is even possible 
to analytically solve the differential equations [121 121] • However, in the present work, 
the flow of the vertex is retained which clearly improves the quality of approximation 

At the end of the fRG flow, the full Green function takes the form = 
iuj5jj> — hj t f(iu) with an effective, noninteracting (but V g -, U- and cu-dependent) 
"Hamiltonian" 

h jd >(iuj) = ho.jj>(itv) - E jV3 v . (19) 

In a last step we have to perform the analytic continuation to the real frequency axis 
iuj — > lo + iO. This is straightforward, as the only frequency dependence of h{iui) is the 
trivial one of the lead contribution Eq. ( |T2l . Then t{oj) can be computed using Eq. ([7]). 

2.3. The NRG approach 

The numerical renormalization group was invented by K. G. Wilson in 1974 as a 
nonperturbative renormalization scheme for the Kondo model [32]. It was later extended 
to the fermionic [231 133] Anderson model which describes a localized electronic state 
coupled to a fermionic bath. The NRG allows thermodynamic and dynamic properties 
of such strongly correlated systems to be calculated at zero and finite temperature 

[22 ESI ESI E3 Ml [26]. 

The key idea of NRG is to discretize the conduction band of the bath 
logarithmically, leading to a tight-binding chain for which the hopping matrix elements 
between the successive sites fall off exponentially with A^q, where A NRG > 1 is the 
discretization parameter, typically 1 < Anrg < 3, and n is the site index. This energy 
scale separation ensures that the problem can be solved iteratively by adding one site at 
a time and diagonalizing the enlarged system at each step, thereby resolving successively 
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smaller and smaller energy scales. Thus, by choosing the length N of the chain so large 
that the corresponding energy scale ~ A-nrg * s smaner than all other energies in the 
problem, all relevant energy scales can be resolved and treated properly. Since the 
dimension of the Hilbert space of the chain increases exponentially with the length of 
the chain, a truncation scheme has to be adopted, according to which only the lowest 
A^kept eigenstates of the chain are retained at each iteration. Recently, it was shown 
that by also keeping track of discarded states a complete, but approximate, basis of 
states can be constructed [39] . This can be used to calculate spectral functions which 
rigorously satisfy relevant sum rules [26] . 

In order to obtain the transmission through the dot t{uj) Eq. (J7J) we follow [26] and 
[35] to compute the imaginary part of the local Green functions at temperature T, using 
the Lehmann representation 

e -UJn/T 

lm Q jtj ,(u) = - vr — - — (n\d jtCr \m) (m\dj f)Cr \n) 5(u - [u m - u n )) 

- 7r — - — (n\d) la \m) (m\dj >a \n) 5(uj + [u) m - w n ]), (20) 

with Z = Yl n e ~ 7 'j ^he m any-body eigenstates |n) and eigenenergies LO n . Since these 
are causal functions, the real part can be accessed by performing a Kramers-Kronig 
transformation [40] . 

In the next three sections we present our results. In Sec. [3] for the U — case. In 
Sec. H] we present the generic phase lapse scenario for interacting spin-polarized dots, 
compare to the mean-field results and investigate the role of finite temperatures. Finally, 
in Sec. |5]we study the spinful two-level dot. 



3. Results: Noninteracting dots 



The large number of parameters makes it essential to analyze the transmission for the 
noninteracting case before considering the effect of two-particle interactions. We focus 
on T = 0. A closed expression for and 9(V g ) (for a fixed spin direction) at U — 

can be obtained from Eqs. (jEJ) and (JTj) by replacing £7(0 + iO) by £Jo(0 + iO), 

2 [rf rf 4 + r^ej + 2Vrfrfrfrf £l e 2 l 1/2 



\t(v g )\ 



(rf rf + rfrf - 2s v^rfrfrf - e l£2 ) 2 + {e 1 v 2 + e 2 v 1 



1/2 



arctan 



eira + ear! 



mod7r 



with Tj = Y2i Tj • For later use we define 



r 



3,1 



(21) 



(22) 



(23) 
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For a fixed set of T l - the 5 dependence of |£(V^)| and 9(V g ) is shown in the first columns 
of Figs. [T]and [2] for U/Y = 0.2. The results are qualitatively the same as those obtained 
for U = 0. For generic level-lead couplings ri- the gate voltage dependence of Eq. ( 1211) 
in the limit of small and large 5/Y is dominated by two peaks (of height < 1) and a 
transmission zero. Associated with the transmission zero is a it phase lapse at the same 
gate voltage. The transmission zero (and phase lapse) follows from perfect destructive 
interference at a particular V g and in the limit of a strong asymmetry in the coupling 
of the two levels to the leads, Y 1 <C Y 2 or vice versa, can be understood as resulting 
from a Fano anti- resonance [20]. Apparently the Fano anti-resonance with vanishing 
transmission is robust if one goes away from this limit towards more symmetric level- 
lead couplings. Across each of the transmission resonances 9 increases roughly by ir 
as expected for a Breit-Wigner resonance. Further details of |t(V^)| and 9{V g ) depend 
on s. For s = + the transmission zero (and phase lapse) is located between the two 
conductance peaks for all 5. For 5 — > the resonance peak positions depend on the 
asymmetry of the 1^- and the separation of the peaks is small if the T l j are close to 1-r 
symmetry, that is close to Yj = T^. For 1-r symmetric dots and 5 = the transmission 
zero (and phase lapse) disappears (not shown in the figures). This is an example of a 
submanifold in parameter space with nongeneric behavior. A complete account of such 
cases (which also remain nongeneric for U > 0) is given in Refs. [12] and [13]. As they 
require fine tuning these parameter sets are presumably irrelevant in connection with 
the experiments and we will here only briefly mention results obtained in such cases. 

For s = — and fixed 1^- the position of the transmission zeros and phase lapses with 
respect to the CB peaks is different for small or large 5/Y (see Fig. [2]). At small 5/Y it 
is located between the two conductance peaks, whereas for large 5/Y it lies on one of 
the outer sides of these peaks [9] . In the crossover regime between these limiting cases 
the height of one of the peaks decreases, while the other becomes broader and splits up 
into two resonances separated by a minimum with non-vanishing conductance [see Fig. 
[2] (g) and (j)]. The crossover scale 5 C depends on the choice of Y 1 -. For large 5/Y, \t\ has 
three local maxima, although the height of one of the maxima is significantly smaller 
than the height of the other two (not shown in Fig. [2]). For fixed, asymmetric Y l j and 
5 — > the separation of the two conductance peaks for s = — is significantly larger than 
for s = + [compare Figs. [T] (a) and [2] (a)]. 

It is important to note that for small 5/Y essential features of the universal phase 
lapse regime established in the experiments are already found at U = 0: regardless 
of the sign s for generic T'- (that is with the exception of a few cases with increased 
symmetry) two transmission resonances are separated by a transmission zero and it 
phase lapse. At U = the peak separation is too small and the shape of the V g 
dependence of the transmission and phase close to the peaks is qualitatively different 
from those observed experimentally (namely Lorentzian-like for the magnitude of the 
transmission, s-shaped for the phase). As we show next the latter problems do not arise 
for sufficiently large interaction U, which in particular leads to an increased separation 
U + 5 of the transmission peaks. 
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4. Results: Spin-polarized dots 

4-1. The generic phase lapse scenario 

In Ref. [19] it was shown that fRG and NRG results for |t(V^)| and 0(V g ) agree 
quantitatively up to fairly large U. For a generic set of couplings 7 we present fRG 
data for 9{V g ) and |t(V^)| together with the occupancies of the levels rij for different 
U and 5 in Figs. CD (s — +) and [2] (s = — ). Increasing U/Y the separation of the 
transmission peaks in the limit of small and large 5/Y increases and is eventually 
given by U + 5. Even though this charging effect appears to be straightforward it 
is important to note that in particular the groundstate at small 5/Y is highly correlated. 
This becomes explicit from the mapping of the present problem on a generalized single 
impurity Anderson and Kondo model as discussed in Ref. [13] . An indication of strong 
correlation effects are the correlation-induced resonances of the transmission found in 
Ref. [12], which we briefly mention below. With increasing U, even at small 5/Y the 
gate voltage dependence of 9(V g ) across the transmission resonances becomes s-shaped 
and the resonances more Lorentzian-like [see third columns of Figs. [1] and [2]. Obviously, 
for s = + the transmission zero and phase lapse remain between the two transmission 
peaks for all 5 and U (see Fig. [1]). For s = — this only holds for sufficiently small 
level spacings as, similar to the U = case, with increasing 5/Y a crossover sets in 
to a regime in which the transmission zero and phase lapse are no longer between the 
peaks. Analogously to the {7 = case, the crossover scale 5 C depends on the particular 
choice of r^. As can be seen from the second row of Fig. [2] (the CB peaks at large U 
have still almost equal height), with increasing U/Y, 5 C is pushed towards larger values. 
The Coulomb interaction thus stabilizes the parameter regime of universal phase lapses. 
This shows that the effect of the Coulomb interaction leading to universal ir phase lapses 
between separated CB peaks in a two-level dot is rather straightforward: for small 5 /Y 
the phase lapse and transmission zero are already present at U = 0, and the effect of 
finite U is simply that the CB peaks become well separated because of charging effects. 
They also lead to a Lorentzian-like lineshape of the peaks and an s-like variation of 9 
across them. — The present scenario has to be contrasted to the one obtained for N > 2 
levels discussed in Ref. [19]. The generic appearance of iV — 1 transmission zeros and 
phase lapses separating the transmission peaks at small 5/Y and U = is specific to the 
case with N = 2 levels. For N > 2 the number of transmission zeros and phase lapses 
at U = strongly depends on the parameters and the mechanism leading to universal 7r 
phase lapses at sufficiently large U (at small 5/Y) is much more involved [19]. This shows 
that although important insights can be gained from studying two-levels to achieve a 
complete understanding of the phase lapse scenario it is essential to study dots with 
Coulomb interaction and more than two levels |19j . 

The lineshape of \t\ shows characteristic differences in the limits of small and large 
5/Y. In the universal regime at small 5/Y and for sufficiently large U/Y the two 
CB peaks have equal width of order Y (not Yj) and equal height which is consistent 
with the expectation that at each peak the transport occurs through both bare levels 
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Figure 1. Systematic account of the energy scales U/T and S/T that govern the 
gate voltage V g dependence of the magnitude of the transmission \t\ (black), the 
transmission phase 9 (red) and the level occupancies (green and blue) of a spin- 
polarized two-level dot at T = 0. The parameters are 7 = {0.1,0.3,0.4,0.2} and 
s = +. For better visibility ni/2 were shifted by 1. The depicted behavior is the 
generic one and in particular qualitatively independent of the actual choice of 7 (up 
to certain cases of increased symmetry; for examples see the text). The behavior at 
U/T = 0.2 is qualitatively the same as the one at U = 0. The results were obtained 
using the truncated fRG. 
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Figure 2. The same as in Fig. [TJ but for s = — . 



simultaneously. A similar behavior is observed in the experiments. In the mesoscopic 
regime (S/T ^> 1) the width of the j'th peak is given by Tj and the relative height hj 
by rjyr^, independent of the value of U. This effectively noninteracting lineshape can 
be understood from the gate voltage dependence of the effective level positions e^ =0 at 
the end of the fRG flow. When one level is charged the effective level position of the 
other level is pushed downwards by U . Besides this, the gate voltage dependence of the 
£ A=o rema i ns linear, leading to two transmission peaks at gate voltages £f =0 (V g ) = 
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with separation U + 5, but with the same width and height as for U = 0. The hopping 
between the two effective levels generated in the fRG flow is small and can be neglected. 

Apart from the CB peaks, for sufficiently large U/Y the transmission shows 
additional features at small V g /Y [see Figs. [T] (b), (c), (e) and [2] (c)]. These are the 
correlation-induced resonances mentioned in the introduction, which have been found 
to be most pronounced at 5 = and in this case occur for interactions larger than a 
critical U c which depends on the T l j and s [12], [13] . Their appearance indicates that the 
groundstate at small 5/Y is strongly correlated (as mentioned in Ref. [T2] the correlation- 
induced resonances are not captured by a mean-field analysis; see below). Associated 
with the correlation-induced resonances is a sharp increase of 9 [see Figs. [1] (b) and 
(e)]. At large U/Y the correlation-induced resonances are exponentially (in U/Y) sharp 
features that vanish quickly with increasing T (see below), which might be one of the 
reasons why up to now they have not been observed in experiments. The correlation- 
induced resonances are not directly linked to the universal phase lapse scenario. 

For increasing U/Y at fixed 5/Y and decreasing 5/Y at fixed U/Y we observe an 
increased tendency towards population inversion of the nj. We define that a population 
inversion occurs if (1) n>i(V g ) = n2(V g pl ) at a certain V g pl and (2) one nj has positive 
and one negative slope at V g pl so that the filling of one level causes a tendency for the 
other to empty. For large U/Y and small 5/Y it is mainly the more strongly coupled level 
(in Figs. [1] and El this is the level 2 shown in blue) whose population increases across 
both CB peaks while it is depopulated in between. This behavior is reminiscent of the 
one discussed in the model with a broad and several narrow levels j3], where a relation 
between population inversion and phase lapse behavior was proposed. Remarkably, for 
sufficiently large U/Y, we find population inversion even for small asymmetries Y2/Y1 
(which is only 1.5 in the example of Figs. [1] and E])- We emphasize that despite this 
resemblance to the observation of Ref. [3], the N = 2 model is not appropriate to 
establish a general relation between the appearance of population inversions and 7r 
phase lapses at small 5/Y [TO], [6]. While the latter are already present at U = 0, the 
former only develop with increasing U [compare Figs. [H (a), (d) or El (a), (d) to Figs. CD 
(c), (f) or El (c),(f)]. Note that the gate voltage V g FI at which the population inversion 
occurs is generically not identical to the position of the phase lapse and transmission 
zero [see Figs. CD(b), (c), (f) and El (c), (f)] [IB]- However, for 1-r symmetric Y l j Eq. (J9j) 
ensures that if a population inversion occurs its position is identical to the one of the 
phase lapse and transmission zero. 

As can be seen in Figs. [1] (c), (f) and El (c), (f) for small 5/Y and large U/Y the 
nj show a rather strong gate voltage dependence between the CB peaks. Nevertheless 
the total dot occupancy n\ + n<i is only weakly V g dependent and close to 1 within the 
entire CB valley. This is reminiscent of the plateau-like occupancy in the local moment 
regime of the single impurity Anderson model showing the Kondo effect. As discussed 
in Ref. [13] a relation to this model can indeed be established. 

We note in passing that with the exception of the nongeneric case of 1-r symmetric 
Y l j, s = + and 5 = 0, the rij are continuous functions of V g . 
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4-2. Comparison to mean-field theory 

In Refs. [10] and [UJ Golosov and Gefen (GG) analyze the phase lapse scenario of the 
spin-polarized interacting two-level dot within the mean-field approximation. However, 
they anticipated themselves that quantum fluctuations not captured in the mean-field 
approach could be important. Examples of this had been pointed out already in Refs. 
[5] and [12]. Thus, GG emphasized that the effects of such fluctuations on their results 
need to be studied in subsequent work. The present subsection is devoted to this task. 

GG consider the subspace of level-lead couplings defined by Tf — Tf = rf — . 
Performing a unitary transformation on the dot states the part of the Hamiltonian Eq. 
(JJJ) containing dot operators can be transformed to (see Refs. [JO] and [UJ) 

#dot + -^lead-dot = 2J ^j d j,adj,a + U d\d X dld 2 ~ t(d\d 2 + H.C.) 
3=1,2 



- [c ,l + 44) + c ,r (iJl - 44) + H.c. , (24) 



with the transformed operators and parameters indicated by a hat. The change of 
basis leads to a direct hopping t between the transformed levels. We here focus on the 
relative sign s = —. GG then introduce the two new dimensionless parameters k and a 



as i = — k5/(2V1 — k 2 ) and a = (\ti\ — \i 2 \)/\Jtl + t 2 - Varying k and a GG investigate 
the phase lapse behavior for fixed level spacing (in the new basis) e 2 — E\ = 5 = 0.256r 
and interaction U/T = 6.4. As the parameters in the new basis are rather complicated 
combinations of the original ones, the variation of a and k corresponds to the variation 
of the (within the above specified subspace), 5 and even s. Our simple picture that 
increasing the level spacing 5 of the untransformed model leads from the universal to 
the mesoscopic regime cannot easily made explicit using the parameters of GG. To make 
the comparison of our results to the mean-field study definite we nevertheless follow the 
steps of GG. 

Varying «,k6 [0, 1[ at s = — we move around in the right part of what is called the 
"phase diagram" by GG (Fig. 4 of Ref. [UJ). In Fig. [3] we show the behavior of \t\ and 
9 varying a at fixed n [Figs. [3] (a)-(d)] and k at fixed a [Figs. [3] (e)-(h)], respectively. 
In the first case, upon increasing a at constant k = 0.5, we move from GGs "phase 
2", (red in Fig. 4 of Ref. [UJ), with the transmission zero and phase lapse outside the 
two CB peaks, into "phase 1" (blue in Fig. 4 of Ref. [UJ), with the transmission zero 
and phase lapse between the peaks. The mean-field approximation correctly captures 
the presence of these two regimes. We find a smooth crossover between them (which is 
why we prefer the notion of different "regimes" rather than "phases"): the V g value at 
which the transmission zero and phase lapse occur smoothly crosses from lying outside 
the right CB peak to lying between the two CB peaks. This is similar to the smooth 
crossover we observe in the s = — case of the untransformed model when moving from 
the mesoscopic to the universal regime (see columns of Fig. [2]). 

For fixed a = 0.6 and increasing k we move from "phase 3" (green in Fig. 4 of 
Ref. [UJ) into "phase 1". In contrast to the mean- field approximation where an abrupt 
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Figure 3. Two traces through the mean- field "phase diagram" of GG. Figures (a)-(d) 
are for constant k and different a moving from GGs "phase" 2 into "phase" 1. Figures 
(e)-(h) are for constant a and different k moving from GGs "phase 3" into "phase 
1" . The inset in Fig. (e) shows a zoom in of the gate voltage region around the phase 
lapse. For a detailed comparison to the mean-field results see the text. The results 
were obtained at T — using the truncated fRG. 



transition from "phase 3" to "phase 1" occurs, Figs. [3] (e)-(h) show a rather smooth 
evolution. The mean-field "phase 3" is characterized by discontinuous population 
switching and a phase lapse between the CB peaks which, surprisingly, is smaller 
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than 7r [TO! E]. Furthermore, in this parameter regime the mean- field results show 
no transmission zero, as discussed in Ref. [6]. However, from Figs. [3] (e)-(h) it is 
apparent that "phase 3", is an artifact of the mean-field approximation: upon taking 
fluctuations into account via fRG, the 7r phase lapse and the transmission zero are found 
to remain in the CB valley. The evolution with increasing k is similar to the s = + 
case of the untransformed model when the level spacing is increased at fixed 1^- (see 
the second and third columns of Fig. [1]). "Phase 3" then corresponds to the parameter 
regime with small level spacing and a sizable U /Y in which correlations are of particular 
importance leading e.g. to the correlation-induced resonances [see Figs. [3] (e)-(h)]. That 
the mean-field approximation fails to properly describe this strongly correlated regime 
is not surprising and has been recognized earlier [12] . In Ref. [13] a connection between 
the small S regime of the present model and the local moment (Kondo) regime of the 
single-impurity Anderson model was established. Thus, the artifacts of the mean-field 
approximation for the spinless two-level dot are reincarnations of the the well-known 
artifacts it produces when applied to the Anderson model in the local moment regime. 

Upon taking fluctuations into account, the discontinuities of the n\ii in "phase 3" 
are washed out. The only choice of parameters for which we find discontinuous behavior 
is the one with 1-r symmetric T l j, s = + and 5 = 0, a case which was already identified 
as being nongeneric [12] . Any arbitrarily small deviation from these conditions leads to 
a continuous gate voltage dependence of ni/ 2 - For parameters close to the nongeneric 
point the change in nj at first sight appears to be rather sharp and an extremely high 
resolution in V g is required to identify the behavior as continuous. 

The fact that the mean-field treatment at U > and small level spacings incorrectly 
causes the transmission zero to disappear and the corresponding phase lapse to become 
smaller than tt is its most consequential problem. By generating such features, the mean- 
field treatment masks the very simple scenario that emerges upon properly including 
fluctuations: for increasing U, the 7r phase lapse and transmission zero found for small 
8/Y at U = remain in the CB valley, while the transmission peaks become well 
separated, Lorentzian-like and the phase acquires an s-shape across the resonances. 

4-3. Finite temperatures 

We next investigate how the phase and magnitude of the transmission are affected 
by finite temperatures. To this end, we use the new FDM-NRG algorithm recently 
proposed in Ref. [26]. In Fig. @]we show NRG data for |t(V^)| and 9{V g ) at different T. 
We consider generic level-lead couplings 7 = {0.27,0.33,0.16,0.24}, U/Y — 10, s = ± 
and 8/Y = 0.02 (universal regime) as well as 5/Y = 4 (mesoscopic regime). As expected, 
with increasing temperature the sharp tc phase lapses are gradually smeared out, the 
transmission zero vanishes and the change of 9 at the phase lapse becomes smaller than 
7r. Furthermore, the CB peaks decrease and broaden. For 8/Y = 4 and s = — [Fig. H] 
(d)] the phase lapse lies outside the CB peaks and outside the window of gate voltages 
shown. 
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Figure 4. Temperature dependence of |i(V^)| and 8(V g ) obtained by NRG for 
7 = {0.27,0.33,0.16,0.24}, U/T = 10 for s = ± and in the universal (8/T = 0.02) 
and the mesoscopic (5/T = 4) regime. We used the NRG parameters Anrg = 2.3 and 
N kcpt ~ 512. 



In the mesoscopic regime [Figs. H] (b) and (d)] the explicit temperature dependence 
of the Green function entering Eq. (El) via Eq. (CO) is rather weak and the temperature 
dependence of the CB peaks and the phase lapse can be understood from the behavior 
in the noninteracting model, but with level spacing U + 5. For small T the height hj(T) 
of the j'th CB peak scales as 1 — hj(T)/hj(0) ~ T 2 /T 2 and the width w of the phase 
lapse asw~ T 2 /(5 + U) 2 [9j. The relevant scale for sizable temperature effects in the 
peak height is thus Tj while it is U + 8 in the smearing of the phase lapses. Since we 
have chosen Tj 5 + U a reduction of hj is visible for temperatures at which the phase 
lapse is still fairly sharp [see Fig. @] (b)]. 

Due to the importance of correlation effects at small S/T, the T dependence of 
|t(V^)| and 9{y g ) in the universal regime is different from the noninteracting case. Here 
the explicit temperature dependence of Q is much stronger and cannot be neglected. 
The resulting T dependence of |t ( V^) | and 9{V g ) is shown in Figs. H] (a) and (c). A 
comparison to Fig. H] (b) shows that in the universal regime the smearing of the phase 
lapse sets in at a lower energy scale than in the mesoscopic regime. This scale depends 
on the relative sign s of the level-lead hopping matrix elements [compare Figs. H] (a) 
and (c)]. Furthermore, in contrast to the mesoscopic regime the scales on which the 
CB peaks and the phase lapse are affected by temperature are comparable. A more 
detailed investigation of the temperature dependence in the universal regime, which 
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also discusses the fate of the correlation-induced resonances, is beyond the scope of the 
present work and is left as subject for future studies. 

5. Results: Spin-degenerate dots 

We finally investigate the effect of the spin degree of freedom on the discussed phase 
lapse scenario, at T = 0. In Figs. [5] (s — +, U/Y — 3) and [6] (s = — , U/Y — 4) we 
show fRG data for the evolution of |t(V^)|, 9{V g ) and nj(V g ) (for a fixed spin direction) 
with increasing 5 for a generic 7 = {0.1,0.2,0.5,0.2}. The overall dependence of 9(V g ) 
on 5 is similar to the one observed in the spinless case (compare to Figs. [Hand [2]). In 
particular, the behavior at small 5/Y appears to be almost unaffected by the presence 
of the spin degree of freedom [Figs. [5] (a), (b) and [6] (a), (b)]. For large 5/Y [Figs. [5] 
(d) and [6] (d)] the transmission resonances are located at odd average total filling of the 
two-level dot indicated by shoulders in the nj [Figs. [5] (d) and [6] (d)]. At these fillings 
and for sufficiently large U/Y the Kondo effect is active and the resonances cannot be 
regarded as Lorentzian-like CB peaks. Instead they show a plateau-like shape known 
from the spinful single-level dot (see Ref. [27] and references therein). Across the Kondo 
plateaus of \t\ the s-shaped increase of the phase is interrupted by a shoulder at 9 ps ty/2 
as expected for the Kondo effect [27J. It would be very interesting to study how each 
of the Kondo plateaus of |t| with increasing temperature crosses over to two CB peaks 
and the phase behaves in the generated CB valley. This question is left for future 
investigations. 

The behavior of the phase in the presence of the Kondo effect was experimentally 
investigated at temperatures comparable to the Kondo temperature [21 J, and much 
below the Kondo temperature [22J. As we study the zero temperature case it is proper to 
compare our calculations to the measurements at low temperatures. Indeed Fig. [6] (c) for 
intermediate 5/Y qualitatively reproduces the experimental results at low temperature 
as shown in Fig. 3 c of Ref. [22]. In particular the increase of the phase by more then 
7r and the absence of clearly developed Kondo plateaus are reproduced. 

In Fig. [5] (a) we left out the fRG data around V g = as for these gate voltages 
some of the components of the flowing two-particle vertex become large. This indicates 
the breakdown of our present truncation scheme [241 [25] and the results for \t\, 9 and 
nj become unreliable. For an explicit comparison to NRG data of |i| see Fig. [7J (a). We 
note in passing that for s = + correlation-induced resonances occur also in the model 
with spin [Fig. (a)] [24] . 

In Fig. [7J we compare fRG and NRG results for 1-r symmetric level-lead couplings. 
The computational resources required to obtain NRG data away from 1-r symmetry 
become large and such data are not required for the aim of the present paper. We can 
then use Eq. (jUJ) and must only compute the occupancies nj, which is numerically less 
demanding. For 5 > (as exclusively shown) it is only the absence of the correlation- 
induced resonances for s = + [compare Figs. [5] (a) and [7J (a)] which is different from 
the results of the generic 7 shown in Figs. [5] and El With the exception of the V g ~ 
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Figure 5. Gate voltage V g dependence of \t\ (black), 9 (red) and the level occupancies 
per spin direction (green and blue) of a spinful two-level dot at T = for different S 
obtained by fRG. The parameters are U/T = 3, 7 = {0.1,0.2,0.5,0.2} and s = +. In 
(a), no data are shown around V g = for reasons explained in the text. 
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Figure 6. The same as in Fig. [SJ but for U/T = 4 and s = — 
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Figure 7. Comparison of fRG (solid) and NRG (dashed) data for \t{V g )\, 9{V g ) 
and rijiVg) (per spin direction) in the universal (small S/T) and mesoscopic regime 
(large <5/r) of a spinful two-level dot at T = 0. The parameters are U/T = 3 and 
7 = {0.15,0.15,0.35,0.35}. In (a) the fRG data for 6(V g ) and n,j(V g ) around V g = 
are not shown. Already the unsatisfactory comparison between fRG and NRG for \t\ 
indicate that the fRG becomes unreliable in this regime. For more details on this, see 
the text. For the NRG parameters we used Anrg = 2.5, Ak op t = 1024 for s ~ + and 
A kcpt = 2048 for s = -. 



regime in the case of small 5 and s = + the fRG and NRG data compare quite well. In 
this case we only show the fRG data for \t\ as the results for the phase and occupancies 
become rather erratic. The reason for the breakdown of the currently used truncated 
fRG is explained in Refs. [23] and [25] and is related to the fact that at small S and 
small V g the correlations in effect become extremely large. 

6. Summary 

In the present paper we studied the appearance of phase lapses in an interacting two- 
level quantum dot considering the entire parameter space using NRG and a truncated 
fRG scheme. 

As a starting point we briefly discussed the noninteracting case at temperature 
T = and pointed out that for generic level-lead couplings, that is up to cases with 
increased symmetry, essential features of the universal phase lapse scenario are already 
established at U — 0. For single-particle level spacings 5 small compared to the level 
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broadenings Tj the transmission is characterized by two transmission peaks of equal 
width with a transmission zero and an associated 7r phase lapse between them (universal 
regime). For large 8/T at U = the appearance or not of a transmission zero and 
phase lapse between the two transmission peaks depends on the relative sign s of the 
level- lead couplings (mesoscopic regime). Within a spinless model we have shown that 
the separation of the two transmission peaks increases linearly with the interaction 
U while the 7r phase lapse and transmission zero remain in the valley between them. 
Furthermore, with increasing U the increase of the phase across the peaks takes an 
s-shape and the peaks become Lorentzian-like, thus assuming shapes resembling those 
observed experimentally [Til [T5l [T6] . For s = — and increasing 8, a crossover occurs 
to a regime in which the 7r phase lapse and transmission zero lies outside the two CB 
peaks. The crossover scale 8 C increases with increasing interaction and thus the Coulomb 
repulsion stabilizes the universal phase lapse behavior. We have investigated the relation 
between phase lapses and population inversions of the level occupancies rij. 

We have shown that a mean-field treatment of the present problem correctly 
reproduces certain features of the behavior discussed above, but is not able to produce 
the universal phase lapse scenario at small 8/T due to artifacts of the approximation 
such as a phase lapse by less than 7r [TOf [TTj and a vanishing of the transmission zero 
[6]. Furthermore, the discontinuous gate voltage dependence of the rij found in the 
mean-field approximation turned out to be an artifact. 

Next, we studied how the phase lapse behavior is affect by temperatures T > 0. 
The universal phase lapse (at small 8/T) is smeared out but remains visible for not too 
large T. In the mesoscopic regime with 8/T ^> 1 the smearing of the phase lapse and the 
decrease of the CB peaks can be understood in detail in analogy to the noninteracting 
case. For 8/T <C 1 correlations are more important and a detailed understanding of the 
temperature dependence of the transmission t(V g ) requires further studies. 

The phase lapse behavior in both the universal and mesoscopic regimes is also 
stable if the spin degree of freedom is included. For sufficiently large U/T in this case 
the Kondo effect is active at odd average dot filling, leading to minor modifications of 
the scenario discussed above. In particular, at large 8/T the T = transmission peaks 
are Kondo plateaus rather than Lorentzian-like CB peaks. Across these Kondo plateaus 
the phase shows a shoulder at 6 ~ tt/2. In contrast, the behavior at small 8/T appears 
to be almost unaffected by the spin degree. A study of the combined effect of finite 
temperature and spin is left for future work. 
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